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Abstract — In this paper, an efficient numerical solution for the 
scattering problem of inhomogeneous dielectric rough surfaces is 
presented. The inhomogeneous dielectric random surface repre- 
sents a bare soil surface and is considered to be comprised of 
a large number of randomly positioned dielectric humps of dif- 
ferent sizes, shapes, and dielectric constants above an impedance 
surface. Clods with nonuniform moisture content and rocks are 
modeled by inhomogeneous dielectric bumps and the underlying 
smooth wet soil surface is modeled by an impedance surface. In 
this technique, an efficient numerical solution for the constituent 
dielectric humps over an impedance surface is obtained using 
Green’s function derived by the exact image theory in conjunction 
with the method of moments. The scattered field from a sample 
of the rough surface is obtained by summing the scattered fields 
from all the individual humps of the surface coherently ignoring 
the effect of multiple scattering between the humps. The statistical 
behavior of the scattering coefficient <r° is obtained from the 
calculation of scattered fields of many different realizations of 
the surface. Numerical results are presented for several different 
roughnesses and dielectric constants of the random surfaces. 
The numerical technique is verified by comparing the numerical 
solution with the solution based on the small perturbation method 
and the physical optics model for homogeneous rough surfaces. 
This technique can be used to study the behavior of scattering 
coefficient and phase difference statistics of rough soil surfaces 
for which no analytical solution exists. 

I. Introduction 

I NVESTIGATION of the radar scattering response of natural 
surfaces is an important problem in remote sensing because 
of its potential in retrieving desired physical parameters of the 
surface, namely its soil moisture content and surface rough- 
ness. Soil moisture is a key ingredient of the biochemical cycle 
and an important variable in hydrology and land processes. 
Although the problem of electromagnetic wave scattering 
from random surfaces has been investigated for many years, 
because of its complexity, theoretical solutions exist only 
for simple limiting cases. Among the existing theoretical 
models the small perturbation method (SPM) [1] and the 
Kirchhoff approximation (KA) [2]. can be mentioned, which 
are applicable for homogeneous surfaces over a restricted 
regions of validity. Numerous techniques based on the basic 
assumptions of the SPM and KA have been developed in 
past in an attempt to improve the region of validity of these 
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models; however, they all have the basic limitations of the 
original models [3]. Other theoretical models such as the phase 
perturbation technique [4] and the integral equation method 
[5] are not applicable, for inhomogeneous surfaces and their 
regions of validity have not been fully determined yet. A 
scattering formulation for rough surfaces with variable surface 
impedance exist [6]; however, the solution is not applicable to 
soil surfaces with low dielectric constant. Several numerical 
solutions of the scattering problem have been proposed to 
identify the region of validity and accuracy of these theoretical 
models. Scattering solution for a perfectly conducting random 
surface using the method of moments has been suggested by 
Axline and Fung [7]. In this method a tapered incident field is 
used as the excitation to eliminate the edge effect contribution. 

A numerical solution for homogeneous dielectric random 
surfaces has been reported [8] where again a tapered illumi- 
nation is used to limit the size of the scatterer. The accuracy 
of the numerical solution with tapered illumination decreases 
with increasing incidence angle. To our knowledge, numerical 
scattering solution for inhomogeneous rough surfaces does not 
exist. 

Analysis of microwave backscatter observations by Oh et 
al. [9] reveals that the existing theoretical models cannot 
adequately explain the scattering behavior of soil surfaces. The 
deviation between theoretical predictions and experimental 
data is attributed to three factors. First, the roughness param- 
eters are often outside the region of validity of the theoretical 
models. Second, the autocorrelation functions measured from 
natural surfaces are very complicated and are not Gaussian or 
exponential correlation functions. Finally, the most important 
reason is that the natural surfaces are not homogeneous in most 
cases, i.e., the moisture content is not uniform in depth. The 
top rough layer which includes clods and rocks is usually dry 
and the underlying soil layer is moist and smooth. 

In this paper, we model a soil surface as an inhomoge- 
neous dielectric random surface comprised of a large number 
of randomly positioned two-dimensional dielectric humps of 
different sizes, shapes, and dielectric constants all lying above 
an impedance surface. At microwave frequencies, the moist 
and smooth underlying soil layer can be modeled as an 
impedance surface, and the irregularities above it can be 
treated as dielectric humps of different dielectric constants and 
shapes. For the field scattered by a single dielectric hump 
over an impedance surface, we have an available efficient 
numerical solution that uses the exact image theory for the 


0196-2892/96S05.00 © 1996 IEEE 



426 


TFFF TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING. VOL M. NO 2. MARCH 10% 


Green's function in conjunction with the method of moments 
[10]. In the solution of a single hump, it has been shown 
that the bistatic scattered field is very weak at points in 
close proximity to the impedance surface; thus, the effects 
of multiple scattering between humps can be ignored. In 
this case, the scattered field from a collection of randomly 
positioned dielectric humps can easily be obtained by summing 
the scattered field of all the constituent humps coherently. 
The scattering coefficients (cr°) is obtained by a Monte Carlo 
simulation. 

In Section n, we summarize the procedure for the numer- 
ical solution of a single hump above an impedance surface. 
Section III outlines the procedures used for generating the ran- 
dom surfaces and for evaluating the statistics of the scattered 
field. Numerical results and their comparison with theoretical 
models are presented in Section IV. 

II. Scattering from Individual Humps 

In this section, we briefly review the procedure for the 
numerical solution of scattering from a two-dimensional di- 
electric object above a uniform impedance surface [10]. The 
radiated field for a dipole source above a dissipative half- 
space medium (Green’s function) is usually evaluated using 
the Sommerfeld integral [11]. This infinite integral, in general, 
is highly oscillatory and computationally rather inefficient. In 
[10], the Green’s function of an impedance surface was derived 
in terms of rapidly converging integrals using appropriate in- 
tegral transforms. The scattering problem was then formulated 
by integral equations which were solved numerically using 
the method of moments. 

Suppose a dielectric object, possibly inhomogeneous, is 
located above an impedance surface and is illuminated by 
a plane wave. The incident field E? induces conduction and 
displacement currents in the dielectric object which together 
are known as the polarization current J e . The polarization 
current can be represented in terms of the total electric field 
inside the dielectric object, which is comprised of the incident, 
reflected, and scattered fields denoted by E*, E r , and E 8 , 
respectively. Thus 

J e (p) = -ik 0 Y 0 [e(p) - 1][JB*(?) + ET{p) + E 8 (p)} (1) 

where k 0 — Y 0 = vWmo. and e(p) is the relative 

dielectric constant of the object at the point ~p — xx + yy. The 
fields E i , ET, and E s are, respectively, given by 

E i (p) = {E^hi + El v { ) exp [ik 0 ki -. p] (2) 

ET (p) = ( RhElhr + RvElvr) exp [ik 0 k r • p] (3) 

S s (?)=ifcoZo / (4) 

where E l h and E l v are the horizontal (E-polarized) and vertical 
(H-polarized) components of the incident field, respectively. 
Rh and Ry are die horizontal and vertical Fresnel reflection 
coefficients and G(p , p') is the dyadic Green’s function of the 
problem. 

There is no known exact solution for the integral equation 
given by (1). Hence, an approximate numerical solution of 
this equation must be obtained using the method of moments. 


This is done by dividing the cross section of the dielectric 
structure into A’ c sufficiently small rectangular cells such that 
the dielectric constant and the polarization current over each 
cell can be approximated by constant values. Using the point- 
matching technique, the integral equation can be cast into a 
matrix equation of the following form: 

l-Zxx] [-2 ly] 0 \X x ) [Vr] 

[■Zvx] ! Zyy] 0 [I v ] = [V,] (5) 

o o [z„jj [pyj ml 

where [Z pq ] is the impedance matrix, [l p ] is the unknown 
vector whose entries are the values of the polarization current 
at the center of each cell, and [V p ] is the excitation vector with 
p = x, y, or z. The entries of [V p ] are simply given by 

U Pi n = 2 /«) *“■ l] 

' 2/n) + ^p{ x ni 2/n )] } * P (6) 

and the entries of [ Z pq ] can be evaluated from 


Zpqmn — ^pq^mn “f* Vn) 1] 

/ I In! 2/m) 

JAS m 

+ Grpq( x m 2/n 5 x mi 2/m)] d$m 


where 8 pq and 6 mn are the Kronecker delta functions, and 
p,q = x, y, or 2 . 

Explicit expressions for the elements of the impedance 
matrix are given in [10] where off-diagonal elements are 
obtained by approximating the Green’s function via its Taylor 
series expansion around the midpoint of each cell and then the 
integration over the cell surface is performed analytically. For 
diagonal elements the free-space Green’s function is approx- 
imated by its small argument expansion and then integration 
is performed analytically. 


m. Monte Carlo Simulation 
of Rough Surface Scattering 

The Monte Carlo scattering simulation of rough surfaces 
made up of a finite collection of dielectric humps involves in 
five major steps as shown in Fig. 1. The first step is to choose 
the type (size, shape, and dielectric constant, etc.) and number 
of constituent humps. The second step deals with generating a 
surface sample by positioning a large number of humps with 
a prescribed probability. The third step in this algorithm is 
to compute the inverse impedance matrices for all constituent 
humps using the numerical method explained in the previous 
section. Next the scattered field from the surface is computed 
by coherent summation of the scattered fields from all of the 
humps in the surface sample. Finally the scattering coefficient 
<r° is obtained by repeating the fourth step for a large number 
of independently generated surface samples. 

The types of constituent humps which in addition to their 
probability of occurrence fully characterize the statistics of 
a random surface. Fig. 2 shows the geometry and dielectric 
profiles of different types of dielectric humps which can be 
handled by this algorithm. For example. Fig. 2(a) shows a 
typical hump for a dry clod above a moist and smooth under- 
lying soil layer (€0 < ei < £ 2 ) and Pig- 2(b) shows the same 
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Fig. 1. Flow chart of the Monte Carlo simulation for the rough surface 
problem. 


hump when the clod and underlying layer are both moist (a 
homogeneous surface). The hump itself may be considered to 
be inhomogeneous as shown in Fig. 2(c). Isolated irregularities 
such as rocks above a flat surface can be represented by the 
hump shown in Fig. 2(d) where the bump occupies only a part 
of the total width allocated for an individual hump. When the 
surface is rough with a short correlation length, the geometry 
of the humps are more complicated and two examples are 
shown in Fig. 2(e) and (f). The profiles of Fig. 2(aMe) used 
in this paper are given by the following functional forms 


y(x) = 

W 2 / TtX \ 

W 

W 


~~A COS \w)’ 

~Y~ X ~ 

T’ 



for (aHc) 



(8) 

y{x) - 

-(-#)■ 

—B < x < B, 

B <W, 



for (d) 



(9) 



numerically by the procedure outlined in [8]. In this procedure 
the hump profile is obtained from a sequence of independent 
Gaussian deviates with zero mean and unit variance which 
are correlated by a set of weighting factors derived from the 
desired correlation function. 

Suppose the set of individual humps includes K different 
humps (including size, shape, and dielectric constant) and the 
profiles of the humps in the set are represented by 
i — 1, • • • , K. Then a sequence of random numbers ranging 
from 1 to K, which is generated by a random number 
generator with the prescribed probability distribution, is used 
to position a large number of humps randomly to construct 
a surface sample. If the total number of humps (M) in the 
surface sample is much larger than the number of constituent 
humps (K) and the random number generator has a uniform 
distribution, the probability of the occurrence of each hump 
in the surface will be about M/K. A functional form of the 
generated surface profile can be represented by 

M / m— 1 \ 

»(*) = H /i» (*- Y , w * ) an 

m=l \ |=i J 


and 


y(x)=AF l (x) + BF 2 (x) 


with 


*(«)-«*• (2 -f) (i) m 

F 2 (z) = cos” (H - |) (l-i)’ 
0 < x < W, for (e) 


where i m ,ii £ {1, •••, K} and W it represents the width 
of the hump of the z/th type. The roughness parameters, 
rms height s, correlation length l, and rms slope m, can be 
computed either numerically or analytically from the surface 
profile given in (11). However, the analytical computation is 
only possible for simple functional forms and simple probabil- 
ity distribution functions. Assuming perfect randomness, the 
(10) average height of the surface can be computed from 


where A and B are constants, n and m are integers, and W is 
width of a hump. The set of constituent humps for a surface can 
be constructed by choosing a finite number of parameters in 
the desired functionals and the desired dielectric constants. The 
profile of Fig. 2(f) is very complicated and should be obtained 


y(x) 


i K r w i 

-i£» / 


fi(x)dx 


(12) 


where L = ^2iLi PiW { and p, is the probability of occurrence 
of the hump of zth type. The rms height s and the rms slope 
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m, respectively, can be evaluated from 


s - ([y{x) - 2 /(*)] 2 ) 1/2 


= j jYl J 0 \p*fi(x) - y(x )] 2 dx J 


1/2 


and 


m = 


' dy(x) _ / dy(x) Y 
dx \ dx / 


2 \ »/* 


{ i K r w ' 

l'51 w ' l [^( z )] 2rfx 


1/2 


(13) 


(14) 


Assuming the surface has a Gaussian correlation function, the 
correlation length l in terms of rms height and rms slope is 
given by 


; = ^ s(nrohdghi) 
m( rms slope) 

It is often required to generate a random surface of given 
rms height s and correlation length l. In that case the required 
surface can be obtained by an iterative process where some 
initial values for the hump parameters are chosen. Then 
the roughness parameters are calculated and compared with 
the desired ones. Depending on the difference between the 
calculated s and l and the desired s and l the hump parameters 
are modified and this process is repeated until the difference 
is below a tolerable error. 

Once the set of individual humps for a random surface 
with given s and l is formed the impedance matrices, [Z pq \i, 
i = 1 , • • • , K, can be computed using the method of moments 
described in the previous section. Since the scattered field 
of a hump near the impedance surface is very weak the 
effect of multiple interaction between humps in a surface 
sample can be ignored. Therefore, by inverting and storing 
the impedance matrices of the constituent humps the scattered 
field of any surface sample comprised of M humps (M > 
K ) and for any incidence and observation directions can be 
computed very efficiently. For a given direction of incidence 
the polarization current in the jth hump for the vertical and 
horizontal polarization, respectively, are given by 




Scatter Angie (Degrees) 

(a) 



(b) 

Fig. 3. Bistatic echo width of a squarod-cosine hump of — 15 + »3, 
W = 0.72A, H = 0.07A over an impedance surface of tj = 0.254 — iO.025 
at (a) 8 t ss 0° and (b) = 45° at / = 5 GHz for E- and H-polarizations. 


Ira 

and 

[I,] j=lZ„]T 1 {V']j (17) 

where j £ { 1, • • • , M} and ij £ {1, • • • , K) representing the 
hump of the ith type. The excitation vector [V]j is computed 
from (6) where the position vector p is specified by the 
discretization procedure and the profile function (11). The 
electric polarization current induced inside the surface sample 
can be represented by 

P = x,y,z 


where [1 p]i is the p-polarized current inside the ith hump. The 
radiated far field can be evaluated from 

E S pp = e i(kof, -* /4) Spp, pp=hhor vv ( 19 ) 

where Spp is the far field amplitude given by 

Sftft = 52 Mx„, y n ) Ax n A y n e -‘*° si ” *■*” 

n=l 

. [e" iA:o cos e * Vn + R B (6s)e iko cos 6tVn } (20) 

S vv = ^Y i Ax n Ay„ e- ik » sto *>*” 

n=l 

* Vn) COS 0 S 

. ([e~* fc ° cos e ‘ Vn - R H {9 3 )e iko cos BtVn ) 



\Z X * 1 


— I 

\m] 

i 

[Zyx\ 

\^yy\ . 

1. 

.[V,]. 


(16) 


(18) 
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Incidence Angle (Degrees) Incidence An gle (Degrees) 

(a) (b) 

Fig. 4. Multiple scattering effect on the backscatter echo width of a surface segment consisting of three squared-cosine humps with €, = 15 + i3 at 
/ = 5 GHz for both (a) HH-polarization and (b) W-polarization. 


Vn ) sin 6 a 

• (e~ iA ° cos + R H { O s ) e iko cos ***•)}. 


( 21 ) 


Here N t is the total number of cell in the surface sample. 

The statistical behavior of the scattered field are obtained 
from the scattered fields E pp of many independent surface 
samples. For sufficiently large N s surface samples the inco- 
herent scattering coefficient is computed from 


°r, = 


2^ Tp 

KLZ 



pp = hh , vv 



Ns 

2* 

£ e; pJ 



3 - 1 



( 22 ) 


where L av = (1 /N s ) i Lj, and Lj is the total length of 
jth random surface. 


IV. Numerical Results 

To demonstrate the performance of the technique proposed 
in this paper, we shall use it to compute the scattering for 
some sample surfaces and then compare the results with 
those predicted by the available theoretical scattering models, 
when conditions apply. First, we consider a surface with 
homogeneous dielectric humps as shown in Fig. 2(a). The 
functional form of the humps are given by ( 8 ) where the 
parameters A and W are varied to generate the set of the 
constituent humps. Keeping A as a constant controlling the 
height and varying W, a set of similar humps can be generated. 
A random number generator with output i € {1, K} 
selects the parameter W{ = BXi, where B is a constant 
controlling the width of the humps and A is the wavelength. In 
this example, the hump parameters were chosen according to 
Table I and the random number generator was given a uniform 
distribution with K — 10 . Before presenting the statistical 


scattering behavior of the surface, it is useful to demonstrate 
the validity of the assumption regarding the significance of 
the effects of multiple scattering among the humps. Fig. 3(a) 
and <b) show the bistatic echo width of a squared-cosine 
hump with W - 0.72A, H = 0.07A, ei = 15 + z3 above 
a surface with 77 = 0.254 — i0.025 (which corresponds to 
€2 = 15 + iZ) at 5 GHz when the incidence angle 6 { = 0° 
and 8 i = 45°, respectively. It is shown that the bistatic 
echo widths at the large scatter angles (near the surface) are 
very weak which implies that the effect of multiple scattering 
between humps can be ignored. In order to illustrate the 
effect of multiple scattering, a surface segment consisting of 
three squared-cosine humps with ei = 15 + i3 above an 
impedance surface with p — 0.254 - i0.025 was considered 
(see Fig. 4). Dimensions of three humps are, respectively, 
given by: W x = 0.8A, H x - 0.08A; W 2 - 1 . 0 A, H 2 - 0.1A, 
and W 3 = 0.6A, H$ = 0.06A. The backscatter echo widths of 
the surface segment were computed twice. In one case the 
scattered field was computed from the polarization current 
of isolated humps (ignoring the effect of mutual coupling) 
and in the other case the polarization current of the three- 
hump structure was obtained directly from the method of 
moments solution (including the effect of mutual coupling). 
Fig. 4(a) and (b) show that the effect of multiple scattering 
is negligible for both polarizations. As long as the ratio of 
rms height to correlation length of the surface (s/0 is small, 
this approximation provides accurate results. For most natural 
surfaces s/l < 0.3 which satisfies this condition [9]. However, 
if the ratio (s/l) is relatively large, the hump type of Fig. 2(f) 
must be used to include the effect of multiple scattering at the 
expense of computation time. 

The rms surface height s and the rms surface slope m for 
this surface can be computed from (13) and (14), respectively. 
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TABLE I 

Roughness Parameters Correspondinc to Constants .4 and B 


Case 

A 

B 

Approx, f 
s l 

in cm 

Exact | 
8 l 

in cm 

At 5 GHz 
ks kl 

Remarks 

1 

30 

0.20 

0.115 2.21 

0.115 2.03 

0.12 2.13 

SPM region 

2 

30 

0.36 

0.208 3.98 

0.207 3.63 

0.22 3.80 


3 

30 

0.70 

0.405 7.74 

0.405 7.15 

0.42 7.49 

PO region 


f Approximation by equations (23), (24) and (15). 
t Numerical evaluation with 4000 humps, 
s : rms surface height, 
l : correlation length. 


TABLE II 

Constants Used in the Numerical Computations 



Individua 

il hump size 

No. of humps 

Length of 

No. of 

Case 

Wi 

dth 

Heij 

ght 

for each 

surface 

segments 

No. 

min. 

max. 

min. 

max. 

surface 

segment 

for a 


(*) 

w 

(A) 

(A) 

segment 

(A) 

surface 

1 

0.2' 

2.0 

0.0066 

0.066 

40 

44 

100 

2 

0.36 

3.6 

0.012 

0.12 

40 

79 

100 

3 

0.7 

7.0 

0.023 

0.233 

40 

154 

100 


and are given by 


s = 


m = 


1 /3 Wf Wf_ 

Lw (.8 A 2 A y + WiV ) 

L i=l ' ' 


1/2 


7 r 


V2A 


where 


and 


1 K 

= — — y' w< 

2AL W ^ * 


*=i 


K 


Lw = ^2 Wi. 


i=l 


(23) 

(24) 


It should be noted that the rms surface slope m of this surface 
depends only on the constant A. Therefore, for a fixed value 


of A, both the rms height and the correlation length increase at 
the same rate with increasing B. Table I shows several values 
of roughness parameters, s and l, corresponding to different 
values of A and B. 

A random number generator was used to select and position 
4000 squared-cosine humps over the impedance surface (77 = 
0.254 — *0.025). Then this surface was divided into 100 
segments to obtain 100 independent surface samples each 
having 40 humps. The length of the surface segment was 
chosen to be about 44A to 154A depending on the correlation 
length of the surface, corresponding to the size of individual 
humps. Table II summarizes the characteristics of the surfaces 
and their constituent humps used in the examples considered 
in this study. 

The backscattering coefficients for the surface at 5 GHz 
with ks = 0.12 and kl = 2.13 (Case 1 in Tables I and II) 
are computed by the Monte Carlo simulation technique for 
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(a) 



Incidence Angle (Degrees) 

(b) 

Fig. 5. Backscattering coefficient a° of the random surface with 
ks = 0.12, kl — 2.13, and e\ — €2 — 15 + *3 as computed by the SPM 
and the numerical technique, (a) HH-polarization. (b) VV -polarization . 

a homogeneous surface with — €2 = 15 + i3 (Fig. 2(b)), 
and compared with the analytical results based on the SPM. 
The comparisons are shown in Fig. 5(a) and (b). For the SPM 
solution, the scattering coefficient <x° is proportional to the 
roughness spectrum (Fourier transform of the correlation func- 
tion). Both the actual and Gaussian correlation functions are 
used in the calculation of the backscattering coefficients using 
the SPM. It is shown that the Monte Carlo simulation agrees 
very well with the SPM prediction when the actual correlation 
function is used. The discrepancies between the Monte Carlo 
simulation and the SPM with Gaussian correlation function 
indicate the importance of the tail section of the correlation 
function in the estimation of < 7 °. 

The numerical simulation was also performed for a surface 
at 5 GHz with ks — 0.42, kl = 7.49 (Case 3 in Tables I 
and II), and €1 = 62 = 15 + i3. The roughness parameters 
of this surface fall within the validity region of the physical 



Incidence Angle (Degrees) 

Fig. 6. Backscattering coefficient <r° of the random surface with 
ks — 0.42, kl — 7.49, and €1 = 62 — 15 + «3 as computed by the PO 
model and the numerical technique for HH-polarization. 
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Fig. 7. The sensitivity of the backscattering coefficient cr° to the dielectric 
constant in case of ks — 0.22, kl = 3.8, and €2 — 15 + t'3 at & = 44°. 

optics (PO) model; therefore, the numerical solution can be 
compared with the PO solution. The scattering coefficient a% h 
predicted by the PO model using the actual correlation function 
agrees very well with the results computed by the numerical 
technique (see Fig. 6 ). In this figure the PO solution using a 
Gaussian correlation function with the same correlation length 
as the actual correlation function is also compared with the 
numerical simulation. It is shown that the agreement is good 
only for low incidence angles ( 0 t < 20 °) and the discrepancy 
between the two solutions becomes rather significant for higher 
incidence angles. In this case, similar to the previous case 
(SPM), it is shown that the tail of the correlation function 
plays an important role in determining the angular patterns of 
the backscattering coefficients. 

With the success of the Monte Carlo simulation in predicting 
the scattering behavior of rough surfaces in the small perturba- 
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tion and physical optics regions, the numerical model can be 
used to study complex surfaces with intermediate roughness 
parameters and inhomogeneous dielectric profiles. For exam- 
ple. the sensitivity of radar backscatter to the moisture content 
of the top layer for the Case 2 random surface is shown in 
Fig. 7 at 0 — 44°. 

V. Conclusion 

In this paper an efficient Monte Carlo simulation technique 
is proposed for computing electromagnetic scattering by inho- 
mogeneous one-dimensional rough surfaces. The surface irreg- 
ularities are represented by inhomogeneous dielectric humps 
of different shapes and the underlying layer is represented 
by an impedance surface. A moment-method procedure, in 
conjunction with the exact image theory, is used for calculation 
of the field scattered by the dielectric humps. It was shown 
that the scattered field near the impedance surface is weak, 
and hence, the effect of multiple scattering between humps 
can be ignored. 

To check the validity of the Monte Carlo simulation, the 
numerical results were compared with the existing analytical 
solutions for surfaces at extreme roughness conditions. A 
smooth surface that satisfies the validity region of the SPM 
and a surface that satisfies the validity region of the PO model 
were considered, and in both cases excellent agreements were 
obtained between the analytical results and those computed 
using the proposed technique. It was found that away from 
normal incidence, the tail of the correlation function plays 
an important role in die determination of the backscattering 
coefficients. 

The analysis presented in this paper is only for one- 
dimensional surfaces and therefore is incapable of predicting 
the cross-polarized scattering coefficients. A numerical 
simulation for a two-dimensional rough surface using a similar 
method is computationally tractable. 
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